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Inclusive production of TJ/if pair in proton-proton interation at LHCb is con¬ 
sidered. This process is forbidden at leading order of perturbation theory, so such 
channels as double parton scattering, XbXc pair production with subsequent radiative 
decays of P- wave quarkonia, contributions of color-octet states, and NLO corrections 
are studied in details. For all these channels we present theoretical predictions of 
total cross sections at LHCb and distributions over different kinematical variables. 
According to presented in the paper results, double parton interaction gives main 
contribution to the cross section of the considered reaction. 
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I. INTRODUCTION 


Inclusive double production of heavy quarkonia was observed in NA3 experiment, where 
production of two J/if} mesons in ttN interaction was studied. The cross section of double 
J/'ip production was suppressed by about three orders of magnitude in comparison with single 
JM. Such suppression can be easily explained at leading order (LO) of perturbative QCD 
lj, 2], In low energy hadronic collisions quark-antiquark subprocess qq —* J/'ifJ/'if gives 
main contribution to the considered reaction, while in the case of high energy interactions 
gluonic subprocesses gg -A J/^J/if are dominant. Leading order O(cc^) QCD predictions 
and subsequent next to leading order (NLO) corrections [3] are in good agreement with 
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experimental results, obtained at LHC jd-9]. It is clear, on the other hand, that at LHC 
energies the gluoinc luminosity is large, so production of two J/ij) mesons in two independent 
partonic interactions (so called double parton scattering, DPS) is possible. In papers 0 ,Q 
it is shown, that cross sections calculated in DPS approximation are comparable with the 
ones obtained in single parton scattering (SPS) approximation. The drawback of DPS is poor 
knowledge of double parton distribution functions. If one accepts the simplest model and 
assumes, that this function is just the product of two single-parton distribution functions, 
the following simple relation holds: 


AB 

'dps 


a, 


SPS°SPS 

a eff 


B 


Using values of single particle production cross sections cr^po and effective cross section 


a e ff ~ 14.5 nb from CDF and DO data 
double J/il) production via DPS: 


12 


14J, it is easy to calculate the cross section of 


a 


j/pj/p 

DPS 


2 nb. 


This value can be compared with SPS cross section, that is calculated at LO 0(c>d) pertur¬ 
bative QCD: 


< 7 , 


JV’Jp 

SPS 


4 nb. 


As for double T(IS') meson production, presented in paper [15| estimates show, that 
cross section of this process is about three orders of magnitude smaller than cross section of 
double J/ij; production. For energy y/s = 7 TeV at LHCb we have 


4?s ~ 8.7 pb, crops » 0.4 pb. 

It can be clearly seen, that while for double J/ij) production DPS and SPS cross sections are 
comparable, in the case of T(1S)T(1S) final state DPS cross section is suppressed by more 
than order of magnitude. 

Question concerning situation in TJ/'ij) production arises naturally. At leading 0(af) 
order in color-singlet approximation this process is forbidden, so the following relation holds: 




a, 


TJ/p 

SPS 


< aJ P T s < a, 


j/ipj/ip 

SPS 


( 1 ) 


The reason is that at leading Oja'/) order for J/ if) J/'ij) and TT processes diagrams with 
f-channel quark exchange are allowed, that saturate low invariant mass region and lead to 
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FIG. 1: Feynman diagrams of gg —>• XcXbSubprocess 

rapid decrease of the cross section with the increase of this invariant mass, a r \j 1/s 2 . At the 
same order of magnitude the contributions to T J/ip production are absent, but production 
via radiative decays of XbXc state is allowed. So, it is interesting to study the modification 
of relation dTJ) with these processes taken into account. 

In the next section cross sections of XbXc production are considered. Section III is devoted 
to calculation of NLO 0(a®) corrections and color-octet contributions to the cross section 
of the process gg —> TJ/'ijj. In section IV we estimate cross section of inclusive TJ/ip 
production at LHCb and compare them with DPS contribution. 


II. gg ->• XcXb 


In the framework of partonic model the cross section of inclusive double heavy quarkonia 
production Q 1 Q 2 is written as a convolution of partonic distribution functions (at LHC 
energies gloun-gluon interaction gives main contribution) and cross section of the hard sub¬ 
process. The latter can be calculated using NRQCD formalism 16], when the cross section 
is written as a series over a relative quark velocity (in the following we will restrict ourself 
to leading order). In such approach the fock column of the heavy quarkonium contains both 
colour singlet (CS) and colour octet (CO) components. According to NRQCD scaling rules 
contributions of CS and CO mechanisms should be comparable, but according to presented 
in 


17 


20] analysis, where experimental results on single Xc meson production at LHC were 


considered, CS contributions dominate. For this reason in the following we will neglect CO 
components of P -wave heavy quarkonia. 
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At the lading order of perturbation theory hard cross sections of the considered processes 
are described by shown in Fig {1] diagrams and the corresponding amplitude can be written 
in the form 


M (99 Q 1 Q 2 ) = { 7^2° (kuPi - fci) Af a) (k 2} p 2 - k 2 ) + 


(k 2 ,pi - k 2 )A^ 2) (. k h p 2 -h)\ , 


u 


( 2 ) 


where fc 1)2 , p 12 are momenta of initial gluons and final quarkonia respectively, s = (ki + k 2 ) 2 , 
t — (pi — ki) 2 , u = {h — P 2) 2 are Mandelstam variables of the partonic subprocess, ei )2 are 
polarisation vectors of initial gluons, and A^(k,q) is the effective vertex that describe 
quarkonium Q production in fusion of two (possibly virtual) gluons with momenta k, q. 
Using NRQCD formalism it is easy to obtain the expressions for these vertices. In the case 
of scalar and axial quarkonia, for example, we have 


A? A. 9 ) 


'2 ( 0 * Q ) 


9% 


ab 


A%'\k,q) = 


3M 3 (M 2 - k 2 - q 2 f 
{- 4 k 2 q a qP - 4 q 2 k a k^ + 2 ( 3 M 2 - k 2 - q 2 ) (,k a q 0 + k?q a ) 

g a13 (( k 2 - q 2 ) 2 + M 2 (3 M 2 - 4k 2 - 4 q 2 )) , 

4 yUU) g% b 

3M 5 /2 (M 2 - k 2 - q 2 ) 

{ (k 2 - q 2 ) ( M 2 -k 2 - q 2 ) (e a ^ u k u + e a ^ u q v ) - 
2 [( M 2 - k 2 - q 2 ) k a - ( k 2 + q 2 ) q a ] e^k v q^ + 

2 [(M - k 2 - q 2 ) q & - (k 2 + q 2 ) k p ] eA^'Kq,,} , 


( 3 ) 


( 4 ) 


where color-singlet NRQCD matrix element 

UR = ^ I^qWI 2 P) 

was introduced. Its numerical value can be obtained, for example, with the help of potential 
models or from the analysis of available experimental data. In the case of tensor quarkonium 
the relation similar to ([3]), (J4|) was also obtained, but it is too complex to present it here. 
It should be noted, however, that in the case of axial meson the effective vertex tends to 
zero in the limit k 2 = q 2 = 0. The reason for such behaviour is that, in agreement with 
Landau-Yang theorem, interaction of axial meson with two massless gluon is forbidden. 
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From presented above amplitude it is easy to get the cross section of the partonic reaction: 

1 


da 

dt 


1287TS 2 


I Mr- 


These calculations were performed in REDBERRY computer algebra system 21] and the 


resulting expressions are, unfortunately, too large to present them here (in computer readable 
form they can be found in supplementary hies for this arXiv preprint). One can, however, 
make a simple estimate in high energy region y/s MjMr- From relation (J2]) it is clear 
that configurations with small transverse momentum pr ~ \Jtu/s give main contributions 
in this region, so, for example, for small t one can take only Erst term in the amplitude (J2]). 
The squared amplitude in this approximation can be written in the form 


where notation 


\M \ 2 « ^B^‘\k 1 ,p 1 ~k 1 )B^\h,p 2 - k 2 ), 


Bi?(k,q) = 

pol 


( 6 ) 


(7) 


is used. If we use the leading term over 1/t and set t ~ 0 in expressions ([7]) (in high energy 
region this is allowed kinematically), in the case of tensor and axial quarkonia we have 

4 M XQ2) ^ 


=| B% Q2 \k,q) = 


64 (0* Q )gl ( _Wq v + k^ 


3 M 3 


(kq) 


( 8 ) 


while for axial meson B^ X1 ' ~ 0. Using relations ((HJ) it is easy to obtain we. 


for scalar and tensor quarkonia decay widths into two massless gluons 


22 ]: 


T (xqo ->• 2 g) = 96cc 


2 |-^oq(0)| 


M 4 


w 128 2 \R' qq (0)\ 

T (xq 2 -t 2 g) - —a s - 


known formulas 

(9) 
( 10 ) 


On the basis of presented above considerations one can obtain estimates for cross sections 
of scalar and tensor quarkonia production in high energy limit (see also 23]): 

r ( Xcj c —> 2 g) r (xbJi, 2 g) 


M 


Xc 


M 


Xb 


d ( 99 XbJ b Xcj c ) ~ (2 J c + 1) (2 Jb + 1) 

The ratios of these cross sections are equal to 


16 4 4 

d (Xc2Xb2) ■ d {XcoXb2) ■ d (Xc2Xbo) ■ d (XcoX&o) = y : y : y : 


( 11 ) 
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If there is at least one axial meson in the final state, the corresponding cross sections should 
be suppressed. 

For numerical calculations we should determine the value of NRQCD matrix element 
(Of). At leading order over relative quark velocity one can use, for example, the derivative 
of quarkonia wave function at origin (see eq. flSl) ), that can be determined from solution of 
potential model. Alternatively one can use relations ([9]), (TTOl) to extract this parameter from 
experimentally known widths of scalar and tensor \q mesons. For charmonia both methods 


give the same value |i?(, c (0)| 2 ~ 0.075 GeV 5 . Presented in papers 181420] analysis, however, 


shows, that experimental data are better described with a larger value 

|R( c (0)| 2 = 0.35 ±0.05 GeV 5 . 


For bottomonium mesons the same approach gives 


171 


\B' 


bb 


1.0 GeV 5 . 


These numbers will be used in the following. We will also neglect the fine splitting of heavy 
quarkonia masses, the following values will be used: 


M Xc = 3.5 GeV, M Xb = 9.9 GeV. 


The scale of the strong coupling constant ct s (/i 2 ) is chosen to be 


/i 2 = 


M 2 + M 2 

Xc ^ Xb 


Using presented above values it is easy to obtain predictions for partonic cross sections of 
X meson production in gluon-gluon interaction. Integrated over total phase space partonic 
cross sections of the considered reactions at energy 


s — s 0 — 2 (M Xc ± M Xb ) 2 , (12) 

for example, are listed in table |H It is clearly seen that for ratios of scalar and tensor mesons 
production cross sections relations (TITj) are satisfied with pretty good accuracy. Dependence 
of some of these cross sections on \/§ and their distributions over t are shown in Figj2j One 
can see, that in the case of tensor meson production in low t, u regions an increase caused by 
f-channel gluon is observed. For axial meson production, on the other hand, in agreement 
with Landau-Yang theorem such increase is absent. 
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Qi/ 0,2 

XcO 

Xcl 

Xc2 

XbO 

16.3 

14.8 

21.4 

Xbi 

2.1 

4.6 

3.8 

Xb2 

21.4 

19.6 

29.2 


TABLE I: Total cross sections of the partonic reactions gg —>• XcJ c Xbj b (i n fb) at s = sq 


(a) 




(b) 



FIG. 2: Dependence of partonic cross sections on \/I (left panel) and t dependence at 'hats = s o 
(right panel) for gg —>• XbXc reactions. On both figures solid, dashed and dotted lines correspond 
to gg ->• XdXbi, 99 XdXb 2 , and gg ->• Xc 2 Xb 2 subprocesses respectively. 


After radiative decays of x,c,b mesons into vector quarkonia reactions gg —* XbXc give 
contributions into cross sections of TJ/ip production: 

da (gg -A TJ/-0 + 2 y) = ^ Br (xbJ h ->■ T 7 ) Br (y cJc ->■ J/^ 7 ) dd (33 -A XMcJc ), 

Jb,Jc=0,2 

where branching fractions of radiative decays are equal to 24] 


Br (% c0 -A J/V»7) = L3 %, 
Br (xd -A J/'tp'i) = 34%, 
Br (x c2 -A J/V»7) = 19 %, 


Br (xft 0 -A T( 1 S') 7 ) = 1.8%, 
Br (xbi -A T(1S') 7 ) = 34%, 
Br (xb 2 -A T(1S') 7 ) = 19%. 


Since branching fractions of scalar mesons’ radiative decays are small, main contributions 
come from axial and tensor particles. At partonic energies s = so presented above values 
give 


d (gg -> XbXc ->■ TJ/-0 + 27 ) = 3.4 fb. 























(a) 


16 18 y/s^ 20 22 24 26 28 



yj$, GeV 


(b) 



l GeV 2 


FIG. 3: Dependence of cross section of T.//'b production in partonic subprocess (left panel) and 
its t-distribution at s = so (right panel). Solid, dotted, and dashed lines correspond to XbXc NLO, 
and color-octet channels respectively. 





T 


FIG. 4: Typical Feynman diagram for gg —>• T J/'i/’ reaction 


Dependence of gg —>■ x&Xc —1 > TJf'ijj + 2'y reaction on partonic energy and t dependence at 
s = So are shown in Figj3] by solid line. 


III. HIGHER ORDER CORRECTIONS 


In the Introduction it was mentioned already that one can bypass the prohibition of 
two vector quarkonia production in gluon-gluon interaction if higher order corrections are 


considered. Typical Feynman diagram describing the process is shown in 


case of gg —> 2J/-0 reaction in was first studied in paper 23] (see also 


25 


ig]4] and in the 


26]). According 
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to these articles the cross section of the reaction can be written in the form 


dd(gg -»■ TJ/ip) _ 257r 3 a^/|/| r 


dt 


46656 m^mf. L 


F{i) + F(u) 


2 


(13) 


where 

m = 

and notations 


rn b rn c 


7r 


d 2 k 


k 2 (k — q) 2 ( 1 + r, 


T~b,c — 


1 + T c r 2 


l + n 1 + T b r 2 


(14) 


4m? r ’ 

b.c 


r = 


2k — q 

|q| 


are introduced. One can easily see that after substitution m b = m c , e b = e c these results 


agree with those presented in ref. 


26]. The integral (fl4|) was calculated numerically, and we 


found that in the domain 0 < |t| < 10 4 GeV 2 it can with pretty good accuracy approximated 
by the relation 

Ae at 


m 


1 - t/Mf 


where 


A = 0.8, a = 2.9 x 10" 3 GeV" 2 , M 2 = 3.6 GeV 2 . 

The mesonic constants in (fT3|) are equal to 

= 0.4 GeV, f r = 0.66 GeV. 

Using these parameters we have obtained the number 

Vsp L s° (99 = 21 fb. 

for integrated partonic cross section at s = sq. Dependence of this cross section on \/I and 
its t distribution at s = §o are shown in Figj3] by dotted line. It should be noted however, 
that these results were obtained in the limit s » |i| 3> Tf 2 T ,, so they can be treated only 
as rough estimates. In our future works we are going to remove this limitation and obtain 
predictions for shown in FigQ] cross sections in full kinematical region. 

The prohibition for TJ/i/j production in gluonic interaction can be bypassed also if one 
consider production of colour octet (CO) states 


27]. According to presented in this article 
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results, main contributions in this approximation comes from gg —y bb ( 3 S [^) cc(' s S ^), gg —> 
bb( :i Sf^)cc( :i S\ 1 ^), and gg — > 66( 3 S'f')cc( 3 S'{ 8 ') subprocesses. Using presented in [27] analytical 
expressions for the cross sections and numerical values 


(pi^) = 1.3 GeV 3 , {Oj) = 9.2 GeV 3 , 

(o^ = 3.9 x 10" 3 GeV 3 , (Oj) = 0.15 GeV 3 , 


for NRQCD matrix elements 28j-|31], it is easy to obtain distribution shown by dashed line 
in Fig]3l According to this figure \/I dependence of different cross sections is quite different, 
with CO channel giving dominant contribution in region of small invariant mass of the final 
pair. 


IV. HADRONIC CROSS SECTIONS 

Let us now consider inclusive TJ/i/i pair production in proton-proton interaction at LHC 
energies. The cross section of this process can be easily obtained from presented in previous 
sections expressions: 

l 

dasps (pp QbQc + X ) = / dxidx 2 f g (xi; g 2 ) f g (x 2 ; g 2 ) da (gg ->■ Q c Q b ), (15) 


where X\^ are gluons’ momentum fractions, and f g (xi^]U 2 ) are distribution functions of 
these gluons in initial hadrons (parametrisation CTEQ6L 32] will be used in the following), 
calculated at the scale g 2 = (M 2 c + M 2 b )/ 2. In LHCb detector conditions, when restrictions 
2 < y cc ,bb <4.5 are set on the rapidities of final quarkonia, used in our paper model gives the 
following predictions for cross sections of T J/ip pair production at \/s = 8 TeV in different 
channels: 


o’sps (PP -t XbXc + X ->■ YJ/ip + X) = 0.2 pb, 
<xsps (pp -> YJ/'i/j + X , NLO) = 1.5 pb, 
cCps (PP -A YJ/ip, CO) = 11.1 pb. 


As it was already mentioned in the Introduction, at LHC energies double parton scattering 
(DPS) could give noticeable contribution to cross section of the considered process. In this 
mechanism two final quarkonia are produced in two separate partonic reactions, and in the 
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FIG. 5: Distribution of pp —>• TJ/U + X reaction cross section at LHCb over TJ/^-pair invariant 
mass. Thin solid, dashed, dotted, and dash-dotted lines correspond to contributions of XbXc, CO, 
NLO, and DPS channels respectively, while thick solid line stands for total cross section 


simplest approximation it can be assumed, that these reactions do not depend on each other. 
Under these assumptions the cross section of T.7/?/; pair production is expressed through 
the cross sections of separate heavy quarkonia production in surprisingly simple relation: 


c’dps (TJ/^) = 


Q'sps('C)o'sps(^/'0) 

a eff 


Using experimental results cr S p S (T) = 0.14 uU cr SPS (J/'0) = 1.28 fib for SPS cross sections 
and 33] and the value cr e // = 14 mb [III, llJ, 341 ]. one can obtain the result 


Cops (pp T J/ijj + X) = 12.5 pb. 


It is clear, that this simple mechanism gives main contribution to inclusive T J/ip pro¬ 
duction at LHCb. Analysis of different distributions, however, could help us to separate 
contributions of other channels. In Figs]5] [6] theoretical predictions for distributions over 
invariant mass of T J /?/; pair and transverse momentum of final quarkonia are shown. On 
both figures thin solid, dashed, dotted, and dash-dotted lines correspond to contributions 





























12 



FIG. 6: Transverse momentum distribution of pp — > TJ/ip + X reaction cross section at LHCb. 
Notations are same as in Figf5] 


of XbXc, CO, NLO, and DPS channels respectively, while thick solid line stands for total 
cross section. For estimation of DPS distributions we have used Pythia8 generator 35] 
with default settings to simulate inpedendent production of T and J/ip mesons and simply 
combined the resulting events. 


V. CONCLUSION 

Presented article is devoted to theoretical analysis of inclusive production of heavy 
quarkonia pair TJ/'ijj at LHCb. It is clear, that at leading order of perturbation theory 
in colour singlet approximation this process is forbidden, so other mechanisms should be 
studied. In our paper we consider such channels as double parton scattering, next to leading 
order corrections, colour octet model and production of vector quarkonia states via radiative 
decays of Xb,c mesons. For all these processes theoretical predictions of total cross sections 
under LHCb detector conditions and distributions over kinematical variables are presented. 

According to obtained in our article results, double parton scattering gives main contri¬ 
bution to the cross section of the considered process. From analysis of different distributions. 
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however, it is also possible to separate contributions of other channel, so one can, for ex¬ 
ample, determine the value of colour-octet NRQCD matrix elements. So, we think that 
additional theoretical and experimental analysis of inclusive T J/ip pair production at LHC 
is very interesting and actual task. 
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financially supported by RFBR (grant #14-02-00096 A) and grant of SAEC ” Rosatom” and 
Helmholtz Association. 
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